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We investigate experimentally the influence of current flow through two independent quantum 
point contacts to a nearby double quantum dot realized in a GaAs-AlGaAs heterostructure. The 
observed current through the double quantum dot can be explained in terms of coupling to a bosonic 
bath. The temperature of the bath depends on the power generated by the current flow through the 
quantum point contact. We identify the dominant absorption and emission mechanisms in a double 
quantum dot as an interaction with acoustic phonons. The experiment excludes coupling of a double 
quantum dot to shot-noise generated by quantum point contact as the dominant mechanism. 



I. INTRODUCTION 

Electronic transport through semiconductor double 
quantum dots (DQDs) has been intensively explored for 
nearly two decades lii^i^ The interplay between a dou- 
ble quantum dot and its environment was investigated 
in detail in previous workg^i^iiiSiiSiiiS using microwave 
spectroscopy. Irradiating double quantum dots with mi- 
crowaves results in photon assisted tunneling (PAT)»ii 
The integration of a quantum point contact (QPC) in 
the vicinity of a single quantum dot allowed charge detec- 
tion,^^ which was later implemented in double quantum 
dot systems ji^ii^ 

The novel application of a quantum point contact as a 
source of energy to drive inter-dot electronic transitions 
in a double quantum dot was recently realizedj ^^'^^i^"^ 
These experiments were explained in terms of acoustic 
phonon based energy transfer between the QPC and the 
DQD circuits. The combination of a capacitatively cou- 
pled DQD-QPC system with time resolved charge de- 
tection resulted in a frequency-selective detector for mi- 
crowave radiation. It allows to detect single photons 
emitted by the QPC and absorbed by the DQD.^^ 

Understanding the back-action of a charge sensor on 
a DQD is important for future possible applications in 
quantum information processings^ The possible dom- 
inant mechanisms that lead to QPC-induced inter-dot 
electronic transitions include electron scattering with 
photons^^ and acoustic phonons^^ or shot-noise^^'^^ de- 
pending on the parameter regime investigated. 

In this paper we study back-action of the current flow 
through the QPC detector on a serial double dot. The 
double dot is tuned to an asymmetric regime, where one 
dot is strongly coupled to the source lead, whereas the 
second dot is more weakly coupled to the drain lead. 
Two independent QPCs can be simultaneously used for 
driving the transitions in the DQD. We observe a non- 
additive effect of both QPCs accompanied by the satu- 
ration of the current across the double quantum dot for 
large QPC currents. We explain the measured data in 



the framework of interaction of electrons with acoustic 
phonons. We relate the power emitted by the QPC to 
the temperature of the phononic bath. The experiment 
excludes the possibility of shot-noise being the source of 
inter-dot transitions. 

This paper is organized as follows. In Sec. [ill we de- 
scribe the fabrication of the sample, its electrostatic char- 
acterization and functionality. In Sec. IIIII we present a 
detailed description of the working regime of a DQD and 
QPCs, followed by the results of our measurements of 
current through a DQD using one and two QPCs. We 
discuss the possible interaction mechanisms in Sec. lIVI In 
Sec.|V]we introduce a model based on electron-phonon in- 
teraction and in Sec. IVII we interpret the measured data. 
Section IVlII contains the conclusions. 



II. SAMPLE AND CHARACTERIZATION 

The sample shown in Fig. [TJa) is based on 
a GaAs/Alo.aGao.yAs heterostructure with a two- 
dimensional electron gas (2DEG) 34 nm below the sur- 
face. It was fabricated by double layer local oxidation 
with a scanning force microscope (SFM).^^ The 2DEG 
is depleted below the oxide lines written on the GaAs 
surfaceSi [white lines in Fig.[lja)]. A 4 nm titanium film 
was evaporated and patterned by local oxidation to cre- 
ate mutually isolated top gates [indicated by the dashed 
lines in Fig. [I^a)]. 

The confinement potential produced by the top gates 
and the oxide lines is shown in the contour plot in 
Fig. mb). It was calculated assuming a pinned surface 
model2^ using the lithographic sizes of the gates mea- 
sured after the sample was fabricated. It shows an ap- 
proximately circular symmetry for the dots, with the left 
quantum dot being slightly larger than the right one. The 
color scale is in arbitrary units. 

The structure presented in Fig. [ija) consists of three 
electronic circuits. The first one is formed by two quan- 
tum dots connected in series [marked by the grey (red 
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FIG. 1: (Color online) (a) SFM image of the sample. In-plane 
gates are defined by the thick white lines, titanium oxide lines 
are indicated by dashed lines. Top gates are labeled with 
black letters. The QDs are indicated with grey (red online) 
circles. The white (yellow online) arrows indicate the positive 
direction of the current in the QPCs. (b) Calculated elec- 
trostatic potential in the 2DEG generated by the oxide lines 
and the top gates, (c) Modulated current component propor- 
tional to the transconductance measured through QPC2 for 
Vqpc2=0.5 mV. An AC voltage of 1.5 mV with a frequency 
of 34 Hz was applied to the in-plane gate ipgl and detected in 
the QPC2 circuit using lock-in techniques. The honeycomb 
pattern and four pairs of triple points are visible, (d) Simul- 
taneously measured double dot DC current. The source-drain 
voltage was Vsd=60 fiV. 



online) circles] and connected to source and drain. A neg- 
ative DQD current corresponds to electrons moving from 
source to drain. Each of the other two circuits contains 
a quantum point contact [white (yellow online) solid ar- 
rows]. A negative QPC current means electrons traveling 
through the QPC in the direction of the arrows. 

To sum up, our structure consists of two barriers defin- 
ing quantum point contacts, two quantum dots, two bar- 
riers determining the coupling of the quantum dots to 
the source and drain and the barrier that determines the 
coupling between the quantum dots. In total, this gives 
seven degrees of freedom and there are seven independent 
top-gates used to tune these barriers and the quantum 
dots. The top gates {tpgl and tpg2) are used to tune 
the DQD into a suitable regime. The top gates tqcl and 
tqc2 can tune the transmission of QPCl and QPC2, re- 



spectively. The middle top gate tm controls the coupling 
between the two dots allowing to change smoothly from 
the single dot regime (large dot spread over the area cov- 
ered by the two red circles) to a weakly-coupled double 
dot. The gates ts and td are used to tune the coupling 
of the DQD to source and drain. 

The potential on both sides of QPCl (QPC2) can be 
lifted with respect to the measurement ground, creating 
a mutual gating effect between DQD and QPCl (QPC2). 
These in-plane gates {ipgl for QPCl and ipg2 for QPC2) 
control the number of electrons on the DQD. 

Due to the presence of the metallic top gates, the elec- 
trostatic interaction between electrons in the quantum 
dots and the QPCs is weakened by screening compared 
to semiconductor-only quantum circuits.^'' The large dis- 
tance between the QPC and the double dot (lithographic 
distance 450 nm) further reduces the sensitivity of the 
QPC for detecting electrons passing through the DQD. 

Figures [TJc) and (d) demonstrate the operation of 
QPC2 as a charge detectorJ^ For both QPCs, the one- 
dimensional subband spacing is larger than 3.5 mV as 
estimated from finite bias measurements. In order to use 
QPC2 (QPCl) as a charge read-out, its conductance was 
tuned to e'^/h. A constant voltage of 0.5 mV was apphed 
between the source and drain leads of the QPC, and the 
current was measured. An AC voltage of 1.5 mV applied 
at 34 Hz to the opposite in-plane gate ipgl {ipg2) mod- 
ulated the current through QPC2 (QPCl). This mod- 
ulated signal which is proportional to the transconduc- 
tance was detected with lock-in techniques. The mea- 
surements were performed in a dilution refrigerator at a 
base temperature of 70 mK. 

The resulting stability diagram of the DQD detected 
with QPC2 is shown in Fig. [TJc). The boundaries be- 
tween regions of different ground state charge configura- 
tions of the DQD are clearly visible. In this measure- 
ment, tpg2 is used to change the number of electrons 
in the right dot and ipgl to change the number of elec- 
trons in the left dot. A few charge rearrangements in the 
lower half of the honeycomb induced by the metallic top 
gate tpgl are present. In general, we find that the top 
gate sweeps lead to significantly more charge rearrange- 
ments than sweeps of the in-plane gates. Change of the 
ipg2 potential combined with simultaneous charge detec- 
tion would result in strong detuning of the charge sensor. 
This would shift the operating point far away from the 
sensitive regime. It can be avoided by using tpg2 be- 
cause it has a weaker lever arm on the QPC. The thick 
line in the bottom-left corner of the plot corresponds to 
a resonance in QPC2. 

In Fig.^d) the corresponding DC current through the 
DQD is plotted. It was measured simultaneously with the 
QPC signal presented in the previous paragraph. The 
source-drain voltage applied to the DQD is 60 fj,V. Only 
two pairs of triple points are visible. Similar sets of data 
can be obtained using QPCl as the detector. 
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FIG. 2: (Color online) (a) Double dot current as a function 
of the gate voltages Vipgi and Vipgg. It represents the charge 
stability diagram of the DQD. The dashed lines outline the 
edges of the honeycomb cells.— M and N refer to the number 
of electrons in the left and the right dot, respectively. The red 
line denotes the detuning axis, with zero detuning occurring 
at the triple point. The red crosses correspond to the positions 
of the current maxima. The green line depicts the calculated 
charge configuration diagram with the tunneling coupling t = 
50 /xeV. It was obtained assuming the constant interaction 
model with a finite tunnel coupling.™ The data was taken with 
VsD = 100 /iV applied symmetrically across the DQD, lifting 
the chemical potential in the source lead /xl up and lowering 
the energy of the drain lead /xr. The insets schematically 
show the level alignment at points indicated by the arrows, 
(b) The same region as in (a) but with a 1 mV DC bias voltage 
applied across QPCl and 60 /iV applied symmetrically across 
the DQD. The triangle indicates the region with the induced 
DQD current. The insets show schematic diagrams of the 
electrochemical potentials in the left and the right quantum 
dot along the selected honeycomb boundaries, (c) Calculated 
current through a double dot as a function of gate voltages. 
The calculation refers to the measurement presented in panel 
(a), where no bias voltage was applied to the QPC. The details 
of the calculation are discussed in the text, (d) Calculated 
stability diagram corresponding to the situations in (b), where 
a voltage of 1 mV was applied across the QPC2. The dark 
blue color in (d) was used to mark the region with a negative 
current. The comparison between (b) and (d) is presented in 
the text. 



III. EXPERIMENTAL DATA 

In the following, we concentrate on a single pair of 
triple points where the DQD showed moderate cou- 
pling. Figure[2fa) shows the DC DQD current (Idot) for 
100 /xeV source-drain bias applied across the DQD. The 
inter-dot mutual capacitance Cm estimated from finite 
bias measurements and from the stability diagram assum- 
ing the constant interaction model^ is 8.8 aF, whereas the 
total capacitance of the left dot is Ci — 86 aF and for the 
right dot, C2 — 76 aF. Each dot contains approximately 
15 electrons and the charging energies are about 2 meV. 

The thin dashed lines in Fig. indicate the bound- 
aries of the honeycomb pattern and the numbers in brack- 
ets (M,N) denote the charge population of the left and 
the right quantum dot, respectively. Here, the left dot 
is strongly coupled to the source lead, whereas the right 
dot is weakly coupled to the drain reservoir. During the 
measurement, both QPCs were kept at zero bias. 

The detuning marked by the dashed grey (red online) 
line in Fig. [HJa) is obtained from the capacitance modei^ 
and expressed by the equation 5 — Ej^ — such that the 
total energy of the DQD, i^tot = El + E^ remains con- 
stant. The energies El and En are the single-particle 
energies in the left and the right quantum dot, respec- 
tively. Converting the energies to gate voltages gives: 

<^ = {Oiipgl^R — aipgl^L)'^Vipgl + {aipg2,R " Ol ipg 2, l) ^Vipg2 ■ 

The lever arms Oiipgij and Oiipg2,j are the lever arms of the 
in-plane gates ipg2 and ipg2 on the left [j—L) or the right 
{j—R) dot, respectively. The voltages Vipgi and Vipg2 
are the voltages applied to the gates ipgl and ipg2. The 
lever arms are extracted from measurements at finite bias 
and from the charge stability diagram of the DQD. The 
obtained values are: aipgi^R — 0.048, aipgi^L = 0.021, 
Oiipg2,R = 0-03 and onpg2^L — 0.04. We take zero detuning 
to occur at the triple point. According to the definition 
above, detuning is positive (negative) in the upper-left 
(lower right) part of Fig. [D^a). Two representative en- 
ergy diagrams are shown in the insets. 

In Fig. IHb) the DQD current was measured in the 
same parameter range at a QPC2 bias voltage of 1 mV. 
The bias voltage across DQD was set to 60 ^eV, i.e., 
smaller than the bias voltage applied in Fig. [21^ a). Despite 
that, the current is strongly enhanced along the bound- 
aries (M,N)-^(M+1,N) and (M,N-|-1)^(M-|-1,N-|-1), cor- 
responding to adding an electron to the left dot. The 
enhancement of the current along the honeycomb bound- 
aries is induced by driving a current through QPC2. An- 
other visible feature induced by biasing QPC2 is the fi- 
nite DQD current in the triangle-shaped area indicated in 
Fig. [DJb) that is normally forbidden by Coulomb block- 
ade. 

The following measurements were carried out with the 
QPCs tuned to their first conductance plateau. The over- 
all experimental results do not depend on this operation 
point of the QPC. 

To investigate the influence of the QPC currents on the 
DQD in the triangular region, we tuned the levels in the 
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FIG. 3: Current through the DQD along the detuning hne: 
(a) for different QPCl currents and (b) different QPC2 cur- 
rents. Fits are plotted with solid lines (see text). 



dot along the detuning line depicted as the solid red line 
in Fig.[2la). Fig.[3]shows the dot current versus detuning. 
The black data points in Fig. [3]Ja) were taken with zero 
bias applied to the DQD as well as to QPCl and QPC2. 
No measurable current above the noise level is detected. 
When a DC current of 50 nA is driven through QPCl, an 
asymmetric peak with a maximum of about 125 fA along 
the detuning line is observed (blue points in Fig. [3fa)). 
This effect is strongly enhanced if the current through 
QPCl is further increased to 75 nA (green points) and 
100 nA (red points). All traces cross zero at the triple 
point (zero detuning). 

A similar, but significantly more pronounced effect is 
observed if QPC2 is driven, as shown in Fig.[3l^b). More- 
over, for negative detuning a small negative DQD current 
is observed. QPC2 is more sensitive as a charge readout 
and it has a stronger effect on the DQD. Therefore, we 
conclude that QPC2 is more strongly coupled to the DQD 
than QPCl. 

We have chosen one point on the detuning line corre- 
sponding to 5=200 fieV and swept the QPCl and QPC2 
currents. The results of this measurement are shown in 
Fig. Ufa). The black (red) filled squares correspond to 
positive currents through QPCl (QPC2) swept from 
to 200 nA. The empty black (red) circles are the traces 
recorded while the QPCl (QPC2) current was swept from 
to -200 nA. The QPC induced DQD current is a lit- 
tle larger in the case when the QPCs are driven with 
positive current. This polarity dependence is significant 
and we can exclude that it is due to a gating effect. As 
mentioned before, QPC2 is more strongly coupled to the 
DQD than QPCl. When the QPC current is swept in 
a positive direction (filled symbols in Fig. S]) the DQD 
current starts to level off (the inflection points for are 
lying between 100 and 150 nA in the QPC current axis), 
whereas for negative QPC current directions this effect 
is not clearly visible. 
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FIG. 4: (a) Current through the DQD at fixed detuning of 
200 fieV as a function of QPCl and QPC2 current. The open 
and closed symbols correspond to the and "— " sign in 
the corresponding expression, respectively, (b) Temperature 
of the bosonic bath as a function of the QPCl and QPC2 
current calculated from (a) (see text). 



Another unexpected feature is observed on the green 
traces. The filled (empty) green squares correspond to 
the QPCl current being swept from to 200 (-200) nA 
while the QPC2 current is simultaneously swept from 
to 200 nA. In a simple picture, we would expect that 
the effects of QPCl and QPC2 are independent and they 
add up, but the measurement contradicts this expecta- 
tion. Due to the action of both QPCs the DQD current is 
slightly larger than in the case when only QPC2 is used. 
In addition, there is an unexpected polarity dependence 
with a maximum DQD current for QPCl being swept in 
negative and QPC2 in positive direction. The remaining 
blue filled squares (empty circles) in Fig. [H^a) were ob- 
tained by driving a positive (negative) current through 
QPCl and a negative current through QPC2. 
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FIG. 5: (Color online) (a) The DQD current as a function 
of the QPCl and QPC2 currents. The data was taken at a 
fixed detuning S = 200 /leV. (b) Calculated temperature of 
the bosonic bath Tb. 

The polarity effect in the DQD current is also present 
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in Fig. ^a), where the dot current was plotted versus 
QPCl and QPC2 currents at fixed detuning i5=200 /xeV. 
In this measurement, the lack of additivity of the effects 
induced by the QPCs is even more visible. 



IV. DISCUSSION OF POSSIBLE MECHANISMS 

A mechanism which can induce the current flow 
through the double dot along the detuning line is pre- 
sented in Fig. inila). The driving current through QPCl 
or/and QPC2 is thought to lead to an emission of energy, 
which can be absorbed by the electron in the right dot. 
If the provided energy matches the energy difference A, 
the electron can be excited from the right to the left dot. 
If the electron leaves the DQD through the left lead and 
the next electron tunnels into the right dot through the 
right lead, then the cycle closes and there is a measurable 
current flowing through the double dot. 

An additional enhancement of the DQD current along 
the honeycomb boundaries as observed in Fig. Hfb), in- 
duced by driving a current through a QPC, can be ex- 
plained in a similar way. In the situation shown in the 
upper inset the electron trapped in the right quantum 
dot can absorb energy emitted by the QPC2, tunnel into 
the left dot and leave the DQD system via the left lead. 
The cycle closes when the next electron tunnels into the 
right dot from the right lead. This QPC2 induced pro- 
cess gives an additional contribution to the DQD current. 
This effect is more pronounced in the vicinity of the triple 
point where energy difference between the levels E]^ and 
En in the left and the right dots are small. The lower 
inset of Fig. [Ifb) shows the analogous diagram for the 
situation when the level in the right quantum dot lies 
above the Fermi energy of the leads. Again, the QPC2 
induced process causes the electrons to move from the 
right into the left contact. 

The possible mechanisms of the pumping effect are 
coupling to acoustic or optical phonons, plasmons, pho- 
tons, shot-noise or thermopower effect. Scattering with 
optical phonons is strongly suppressed as long as the rele- 
vant energy scales are smaller than the optical phonon en- 
ergy;?® Coupling to plasmons can be ruled out as welli^ 
We can also exclude the shot-noise as a source of the 
energy, because during the experiment both QPCs were 
tuned to their first plateau. Measurements performed at 
0.5Go and 1.5Go (Go = 2e^//i) showed a qualitatively 
and quantitatively similar behavior. This is in contrast 
to previously measured datal^i2i where no DQD current 
was observed in the plateau regions. Coupling to acous- 
tic phonons is the most likely mechanism of inducing the 
current in the DQD. Further below in this paper, we 
discuss the data in the light of phonon coupling and a 
related thermopower effect. 

The questions arising from the data presented above 
are the following: is the strong difference of the peak 
heights on the positive and negative side of the detuning 
(Fig. [3]) due to the asymmetry in the DQD coupling to 




FIG. 6: (Color online) (a) Driving current through the QPC 
induces a current through the DQD. An electron absorbs an 
energy quantum A emitted from the left or the right QPC 
and is excited to the excited state. It leaves the dot via the 
left barrier and another electron can tunnel from the right 
into the right dot. (b) Schematic energy level diagram for the 
DQD. The wave functions of the bonding (^b) and antibond- 
ing ('I'a) states are shown. The electrochemical potential of 
the left (right) lead is /il (mr)- 

the leads? Why do the effects of the QPCs not add up? 
What is the reason for the saturation of the DQD cur- 
rent observed in Fig. SKa)? What is the mechanism of 
the energy transfer from the QPCs to the DQD? Can it 
explain the polarity dependence? In the next section we 
present a model that attempts to answer most of these 
questions. 

V. THE MODEL 

In the following we derive a scenario, which explains 
the pumping effect based on electron-phonon interaction. 
First, wc introduce the two- level system describing the 
DQD. Then, we consider all possible transitions between 
different energy states of the DQD and express them in 
terms of tunneling rates. Subsequently, we derive the 
energy dependence of the tunneling rates. The intra- 
dot transitions are calculated in a framework of electron- 
phonon interaction. Next, we set up a master equation 
and obtain the complete expression for the DQD current 
as a function of detuning and the temperature of the 
phonon bath. 

Close to a pair of triple points, a double quantum dot 
can be regarded as a two-level system, ^'■^^ whose bonding 
(ground) and antibonding (excited) states Eb and Ea 
are separated by an energy A = \^6^ + At'^ as shown in 
Fig. [6jb),^^ where t is tunneling coupling between the 
dots. The corresponding eigenvectors for bonding and 
antibonding states are |5'b) and |^'a)- The components 
of the bonding and antibonding eigenstates in the basis 
of and \<j>n), the wave functions in the left (L) and 
the right (R) dot, are Ci j = (0i|4'j), where i = L,R and 
J=B,A. 

As in Sec. Illll the detuning is defined as S = El — En 
and the total energy is i?tot = -E-l + E^. We assume. 
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that the number of the electrons in the quantum dots 
is fixed and its ground state energy is -Egso = 0. If 
we add one extra electron called an excess electron the 
ground state energy will be i^csi — ^(^'tot — A), while for 
two excess electrons, the ground state energy is i?GS2 = 
\ {Etot + A) + e'^/Cm- The corresponding electrochemical 
potentials are fii = Eqsi and fi2 = Eqs2 - Eqsi- 

In the vicinity of a pair of triple points a double quan- 
tum dot can have one out of four different charge states. 
These different charge configurations are presented in 
Fig. [T] The "empty" state corresponds to a situation 
where there is no excess electron present in a dot and the 
occupation probability of this state is peso- The index 
GSO denotes the zero-electron ground state. In addition, 
one excess electron may occupy the bonding (ground) 
state with probability pcsi or the antibonding (excited) 
state with probability Pexi- The last possible charge 
configuration is when there are two excess electrons in 
a double quantum dot (two-electron ground state) with 
occupation probability pgs2 ■ 

The transitions between these states are determined 
by the tunneling rates Tij and the thermal broadening 
of the Fermi function fij in the leads. The index i=L,R 
denotes the left (L) or the right (R) barrier through which 
the electron tunnels and the index j=0,l,2,3 labels the 
transition (see Fig. [7]). For example, if the dot is in a 
zero-electron ground state (GSO) and the electron tunnels 
in via the right lead, the corresponding rate is Fft^o/R.Oi 
as shown in Fig. [71 For j—0{l) the Fermi function is 
/i,o(i) = l/{exp [(£'b(a) - fiO/kBTp] + 1} where fee is 
the Boltzmann constant and Tp is the temperature of 
the lead. For j=2,3 the expression is analogous but the 
energy is lifted by the mutual charging energy e'^/Cm- 

In order to explain the experimental data presented 
above, we have to take into account that the tunnel- 
ing rates Tij do depend on the electronic wave function. 



They can be expressed as 



cients c^j are the left (i = L) and right (i 



The coeffi- 
= R) energy 



dependent components of the eigenvector of the wave 
function corresponding to the bonding (j=0,3) and an- 
tibonding {j—1,2) states. The amplitudes ji are energy- 
independent parts of the F^j . 

The rates describing the inter-dot processes, that is 
absorption Fabs and emission Fom arc marked in Fig. [T] 
In the following, we assume that the double quantum 
dot is coupled to a bosonic bath in thermal equilib- 
rium described by the Bose-Einstein distribution func- 
tion nB(£',rB) = 1/ [exp (E'/fceTB) - 1]. The tempera- 
ture Tb of this bath is determined by the current of the 
QPC and the base temperature of the cryostat. In case of 
coupling to acoustic phonons, the emission and absorp- 
tion rates can be expressed as (derivation is presented in 
Appendix A): 

Tem/abs = ("^ (A, ^b) + ^ T ^) 9X (A) , 

(1) 

where the index A denotes piezoelectric transver- 
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FIG. 7: (Color online) Energy level diagrams presenting four 
possible charge states of the DQD in the vicinity of a pair 
of triple points. The state GSO corresponds to the situation 
with no excess electron present in the double dot. The states 
GSl and EXl are the ground and excited state of the dou- 
ble dot in the case of one excess electron, respectively. The 
fourth state GS2 is a two-excess electron ground state. The 
occupation probabilities of these states are labeled as paso, 
Pgsi, Pbxi and pgs2. The transitions between these states 
are determined by tunneling rates Fij and Fermi function 
fi,j. The index i — (L,R) correspond to the transition occur- 
ring through the left i — (L) and through the right i = (R) 
lead. The second index j = 0, 1, 2, 3 labels the transition of 
the electron in the DQD. 



sal phonons (pe,T), piezoelectric longitudinal phonons 
(pe,L) or longitudinal deformation potential coupling 
phonons (dp,L). The exponent n\ is 1 for piezo- 
electric phonons and n\=3 for deformation potential 
coupling. The upper (lower) sign stands for emission 
(absorption) of energy quanta. The values of the energy- 
independent coefficients 7a are given in Appendix]^ The 
form- factor (?a(A) is represented as: 
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where tq denotes the radius of a single QD, d is the dis- 
tance between the dots, ca is the speed of sound of A- 
phonons. The complete expression for the double-dot 
geometry factor h\ {Ad/fic\) is given in Appendix 1X1 In 
Eq. ^ the first factor Cl b^l a related to the symmetry 
of the double quantum dot wave function. It suppresses 
transitions for the asymmetric system. The second term 
of Eq. ([2]) refers to the shape of the individual dots. It 
gives a high energy cut-off for phonon wavelengths much 
smaller than the size of a single dot tq. The last term 
of Eq. ^ arises from the separation of the single dots 
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in the double dot system. It suppresses small energy 
absorption for phonon wavelengths much larger than d. 
Only phonons with a wavelength comparable or larger 
than the DQD separation can interact with the electron. 
For large energies this term has oscillatory behavior^". 

To investigate the influence of the QPCs on the dot 
presented in Fig.[3]in terms of rates and occupation prob- 
abilities defined in Fig. [7] we only take into account the 
processes surrounded by the dashed line. This is reason- 
able, due to a large mutual charging energy relative to 
the tunneling coupling. Counting the electrons passing 
through the right barrier leads to the following expression 
of the current through the double dot (the derivation is 
presented in Appendix [B|) : 

Idot =|e| [— PGSirR2/R2 — PEXirR3/R3 

+PGS2rR2 (1 - /r2) +PGS2rR3 (1 - /rs)] • 

The occupation probabilities posi, Pexi and pgs2 are 
functions of F^ ^ and j resulting from a stationary so- 
lution of the master equation (see Appendix IB|) . 

VI. RESULTS AND INTERPRETATION 

We used expression (jS]) to fit the data shown in Fig. [31 
During the fitting procedure, the 8 traces shown in 
Fig. EJa) and (b) and an additional set of 16 traces being 
a combination of /qpci,2 — 0, ±50, ±75, ±100 nA were 
fitted simultaneously. The following fitting parameters 
were shared: amplitudes 7l, 7r and tunneling coupling 
t. The only parameter specific for each trace was the 
temperature of the phonon bath Tb. 

The fit represents the overall shape of the measured 
data very well. The extracted tunneling coupling is 
t = 50 fieV. In Fig. ^a) the grey (green online) solid 
line is the calculated boundary between the honeycomb 
cells assuming a tunneling coupling of 50 /xeV, the black 
(red online) crosses are the maxima position of the DQD 
current peaks and the black dashed line is the bound- 
ary assuming t = 0. Unfortunately, the sample was not 
stable enough to map a stability diagram with a reso- 
lution high enough to determine the tunneling coupling 
directly. However, the t obtained from the fits seems to 
be reasonable and matches with the data presented on 
the stability diagram. 

The temperature of the phonon bath Tb obtained from 
the fits varies from 0.6 K (blue trace in Fig.[31Ja)) to 1.2 K 
(red trace in Fig. [3I^b)). The difference between the tem- 
peratures Tb and the electronic temperature Tf==100 mK 
gives rise to the DQD current. 

The extracted amplitude 7r is of the order of 7.8 MHz 
and 7l is about 0.5 GHz. This is in agreement with our 
previous statement that the right barrier is more opaque 
than the left one. This leads to a deviation from perfect 
antisymmetry of the dot current along the detuning line, 
i.e., to suppressed dot current for negative detuning. 

The tail of the curve at large detuning is mainly deter- 
mined by the amplitudes gC^ a present in Eq. ^ that 



drop to zero like t'^ /S"^ and the Bose-Einstein distribution 
function given in Eq. ([1]). Although, the real distribution 
of the bosonic environment is not necessarily equilibrated 
and could have another form, the qualitative agreement 
with the data docs not depend strongly on the details of 
this distribution. 

Using the parameters obtained in the fits, we cal- 
culated the current as a function of the gate voltages 
using the model containing all four charge states pre- 
sented in Fig. [71 In the absence of the QPC current 
(Tb = Tf^70 mK), the result is shown in Fig. ^c). 
The asymmetry and magnitude of the current along 
the upper dashed line is reproduced very well. This is 
not the case for the lower branch indicating transitions 
(M, N) ^ (M + 1, N). This may be due to a change of 7l 
and 7r, which we assumed to be constant in our model 
but which may change in the experiment. For 1 mV DC 
bias across the QPC2 (Fig. [2jb)), the bosonic tempera- 
ture is around Tb=:0.8 K. The corresponding calculations 
are shown in Fig.[2Id). The red triangle indicates the re- 
gion with the QPC induced current that is in agreement 
with the measured data in (b). 

To investigate the dependence of the phonon tempera- 
ture on the QPC current we calculated Tb for every point 
from Fig. [H^a) using the values obtained in previous fits. 
The results are shown in Fig. [Hb). For small QPC cur- 
rents the error bars are large and no clear dependence is 
visible. For QPC currents above 50 nA the dependence 
is quadratic Tb ~ ^qpc which means that the temper- 
ature of the bosonic bath is proportional to the power 
emitted by the QPC. The reconstructed temperatures 
corresponding to the measurement presented in Fig. [51^ a) 
are plotted in Fig. [^b). For very small DQD current it 
is impossible to estimate the temperature of the phonon 
bath with sufficient accuracy. 

The saturation of the DQD current for large QPC cur- 
rents cannot be attributed to the high occupation proba- 
bility of the antibonding state pexi • We have estimated 
that the pexi value does not exceed a few percent and 
most of the time the dot is occupied by one electron in 
its ground state. The maximum current is determined 
by the right tunneling barrier. 

A series of experiments reporting the observation of 
a DQD current induced by a single and independently 
biased QPC, is described in Refs. [islfl^fTTl . These ex- 
periments were performed in the regime of large DQD 
current, strong bias voltage applied to the QPC and 
large tunneling coupling, which is similar to our situa- 
tion here. The DQD current was related to inelastic re- 
laxation of electrons in partly transmitting ID channels 
of the QPCl^ and qualitatively consistent with an en- 
ergy transfer mechanism based on nonequilibrium acous- 
tic phononsji^iii In contrast to our experimental data, in 
these experiments the DQD current was large when the 
conductance of the QPC was tuned to e^/ft. and strongly 
suppressed in plateau regions. 

Recent research has proven that the absorption of a 
photon can be the dominant processi^i^ in similar sit- 
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uations. However, these time-resolved experiments were 
performed in a different regime, where the dominant tun- 
neling rates are of order of 1 kHz, whereas in our system 
the double quantum dot is much more strongly coupled 
to the leads. Another difference is the presence of a Ti 
top-gate in our structure, that screens the electrostatic 
interaction between the DQD and the QPCs. Due to the 
lower sensitivity of direct dot current measurements com- 
pared to the time-resolved technique it is not possible to 
observe the gap in the DQD current when leVgpcl < S. 
Calculations of the emission and absorption rates in a 
DQD induced by electron-photon interaction show, that 
the effect is irrelevant compared to the emission and ab- 
sorption of phonons discussed here.— 

We have also tested the possibility that the entire dot 
current in the region forbidden by the Coulomb blockade 
is due to a thermopower effect induced by different tem- 
peratures in the source and the drain. We found, that 
it would be only possible if the temperature difference 
between source and drain lead was larger than IK for 
a QPC current of 100 nA, which is one order of magni- 
tude larger than expectedj21 Another argument against a 
thermopower model is that always the drain lead would 
have to be warmer, even if the far QPCl (that couples 
better to the drain lead) was biased. Even so, the ther- 
mopower model did not describe the data as well as the 
emission/absorption model. 



VII. CONCLUSIONS 

We have presented the influence of two independent 
quantum point contacts on a double quantum dot. A 
number of of questions arising during the investigation 
and presented in Sec. IIVI could be answered. 

In the first place, we established the possible dominant 
mechanism of energy transfer between the double quan- 
tum dot and the quantum point contact in the investi- 
gated regime. Driving current through the QPC leads to 
emission of energy that increases the temperature of the 
bosonic environment. We identify these bosons as acous- 
tic phonons. To model the interaction of the phonons 
with the double quantum dot we have assumed that their 
energy distribution is described by Bose-Einstein statis- 
tics. 

Another important point is a non-additive effect of 
both QPC currents. It is understood in terms of the tem- 
perature of the bosonic bath. We find that the DQD cur- 
rent is proportional to the power emitted by the QPCs. 

Next, we interpreted the leveling off of the DQD cur- 
rent as a QPC current is increased. For large QPC powers 
(above 0.1 nW) the temperature of the phonon bath in- 
creases linearly. The observed saturation of the current 
is due to the finite transparency of the tunneling barriers 
and not to the high occupation probabilities pexi- 

The polarity dependence (Fig. 2]) cannot be explained 
within the discussed model and its origin remains to be 
investigated. It would be interesting to further investi- 



gate this effect, for example, by using different geometri- 
cal arrangements. 

Finally, we observed strong deviation from perfect an- 
tisymmetry of the dot current along the detuning line 
when the QPC current is driven. It can be attributed to 
the asymmetry of the source and drain barriers. 

All the measurements were performed with both QPCs 
tuned to their first plateau. Thus we can exclude the 
influence of shot-noise phenomena in the quantum point 
contacts. 



APPENDIX A: ABSORPTION AND EMISSION 
RATES IN A DQD INDUCED BY 
ELECTRON-PHONON INTERACTION 

Generally, the emission and absorption rates can be 
expressed using Fermi's golden rule: 

rcm/abs - y E I |i^c-p (r) |$A> |' ^ (A - ki,^,^) , 

(Al) 

where the sum extends over all wave vectors q. The index 
A denotes the type of acoustic phonons and their coupling 
in GaAs: piezoelectric longitudinal (pe,L), piezoelectric 
transversal (pe,T) and longitudinal, deformation poten- 
tial coupling (dp,L). The phonons have linear dispersion 
relation u>q.x = c\\<i\- and $a are wave functions 
of bonding and antibonding states separated by energy 
A. The interaction hamiltonian i?c-p can be written as 
a sum of piezoelectric interaction iJpo and deformation 
potential coupling H^p: 



i/e-p (r) - ffpe (r) + ffdp (r 



(A2) 



with 



^^po (r) 



lei di4 



1/2 



F(q)e'' 



1/2 



|q|e 



(A3) 



In above equations du is an element of piezoelectric ten- 
sor, Co is the vacuum permittivity, e is the dielectric con- 
stant, N is the number of atoms in the crystal, M is 
the atomic mass and D denotes deformation potential 
coupling constant. The dimensionless function F (q) has 
form: 

|q| ^ki 

where Siki is the Levi-Civita symbol and e^^i is the i- 
component of the eigenvector associated with mode A. 
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Inserting Eq. (jA2p into Eq. (|Aip leads to the following 
expression: 



- cm/abs 



(7peA [.gpe,T (A) + .gpe,L (A)] + TdpA^gdp (A)) 



X ( n(A) + ^ T ^ 



(A4) 



where the upper sign refers to phonon absorption and the 
lower to phonon emission. The constants 7po and 7dp are 
given by: 



7pc ^ 
7dp 



2tt (2eeo)^ ph'^cl 



5 X 10^" meV~^s-\ 



= 7.3 X 10" meV~^s-\ 



where p is a density of the GaAs crystal. The energy 
dependent functions ^pc (A) and ^dp (A) are defined as: 



5pe(A)-E|^/AFA(±q) 



X 



|($B|e±''in*A)|'(5(A-;ia;q,A) 



5dp (A) 



fC'cl 



47rA4 

X |(<i>B|e±^'i'-|*A)r5(A-;ia;qa). 



(A5) 



Assuming negligible overlap between the wave func- 
tions of the two dots and taking a gaussian-shaped single- 
electron wave function, the matrix element is found to be: 

|(<i>B|e±''i'-|'i>A)r = 2c^,BctAe-'''""°'' [1 - cos (qd)] , 

(A6) 



where d is the distance between the dots and rp is the 
radius of a single dot. Inserting (jA6p into Eqs. (jA5p gives: 



(A) = .e-^'i"-'') h. 



For piezoelectric transversal phonons the above expres- 
sion was calculated by averaging the function /po,T (q) 
over all possible transversal directions. The geometry 
factors h\ {rj = Ad/hc\) are given by: 



24 

hpcL = ^ + ''S'?"^ [9?/ {rf - 10) cos?7 

+ {t]'^ + 39?7^ + 90) sin 77 1 , 
32 

V.T (^) = ^ + 16??"^ [v {V^ - + 405) COST? 



- 3 (377" - 62r/2 -I- 135) sinry] , 
^dp,L (??) = 2 - 2?7"^ sirn]. 



(A8) 



Combining Eq. (|A8|1 with Eq. (|A7|) and inserting the 
result into Eq. (jA4p gives a complete expression for the 
absorption and emission rates. 

APPENDIX B: RATE EQUATION 



To relate the DQD current to the tunneling rates we 
write down the rate equation for the occupation of the 
states: 



d 
dt 



PGSl 




PEXl 




PGS2 





Tabs 

A2 



rem B2 

-(Tem + As) Bs 

A3 -{B2 + B3) 





PGSl 




PEXl 




_PGS2 _ 



(Bl) 



with additional condition posi +Pexi +PGS2 = 1- The positive (negative). 



terms Aj and Bj are defined as: 



Aj — ^i,jfi,j 
i=L,R 

Bj = ^ (1 - /ij ) . 



(B2) 



i=L,R 



To find the expression for the current flowing through a 
DQD, we take the right barrier as a current reference. 
It means that, if an electron passes the right barrier to 
the left (right), its contribution to the DQD current is 



Idot — |e| [— PGSirR2/R2 — PEXiTrs/rs 

+PGS2rR2 (1 - /r2) + PGS2rR3 (1 - fm) 



(B3) 



The first (second) term of Eq. (|B3|) corresponds to the 
electrons moving from the bonding (antibonding) state 
to the right lead and the third and fourth term to the 
electrons entering the ground or excited state of the dot 
from the right lead. By inserting a stationary solution of 
Eq. (jBip into Eq. (jB3[) one obtains an expression for the 
steady state DQD current. 
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